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!>• ■ Abstract 

\o , 

A confined quark model study in a couple of chirally anomalous processes is 

>n 

presented in comparison with effective meson Lagrangian approaches of various 

\ kind. The processes considered are ir° — > 77 and 7 — ► 7r + 7r _ 7r° (or the equivalent 

Oh! 

7T7 — ► 7T7r) for which there is a well-known low energy theorem relating the latter 
amplitude with the former one by a very simple algebraic relation in the zero energy 
(or chiral, or soft pion) limit. Our quark model naturally generates the so-called 
contact term in the amplitude for the second process, but with the opposite sign to 
what effective chiral meson models indicate. A reinterpretation of our vector pole 
contribution restores the consistency, however. While the first reaction is observed 
to serve in calibrating various models, it is found difficult, based upon the quality 
of the existing data in the second reaction, to single out the best model, which 
appears indispensable for testing the validity of the above low energy theorem. Thus 
the proposed experiments and their analyses should aim at attaining an optimal 
accuracy. 
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1 Introduction 



A little more than a quarter of a century ago, the chiral (axial) anomaly was found to 
amend the ordinary PCAC (partial conservation of axial current) relation in order to 
account for the physical decay of the neutral pion into two gammas, which otherwise 
would be too small and vanish in the strict chiral (zero meson mass) limit As has 
been frequently quoted, a quantitative comparison of the theoretical (anomalous PCAC) 
prediction with the measured decay rate lead to one of the strongest evidence for the three- 
colour nature endowed to the quarks. More precisely, together with a proportionality 
constant N c the anomalous PCAC gives the n° — > 77 amplitude 

i^ 77 = aN c /3nF n = 0.025 GeV -1 , (1) 

where F w m 93 MeV is the charged pion decay constant. In comparison with the cor- 
responding quantity extracted from experiment: F^ f = (0.0250 ± 0.0003) GeV -1 , one 
obtains N c = 3: the number of colours, see for example ||. 

A couple of years later, a low energy theorem which relates the F„- 77 to the amplitude 
for the 7 — > 37r reaction at zero energy (more precisely the soft pion and soft photon limit) 
which we shall denote as F 3n (0) = F 3n (0, 0,0) was deduced based upon the anomalous 
PCAC and current algebra by Adler et al. ||, Terent'ev and Aviv and Zee 0: 
henceforth we shall call this low energy theorem as the ATA theorem for brevity. The 
occurrence of this 7 — > 3tt was also recognised as due to chiral anomaly, and the theorem 
claims 



ei^(0) = F^/Fl (2) 

or 

fW(0) = -!= = F~ = 9.54 GeV- 3 . (3) 

Note that the value of the chiral anomaly prediction denoted as F^ om in the above equa- 
tion will appear frequently in what follows. Later it was shown that, from the electro- 
magnetically gauged Wess-Zumino-Witten (WZW) term ||0 involving the chiral fields 
(ordinary pseudo-scalar octet mesons), both F^-y and F^ om can be obtained as the tree- 
level amplitudes explicitly satisfying the ATA theorem almost for free |7]] |§ |§ • 
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To test the ATA theorem, several physical processes like ir —>■ 2tt in the Coulomb field 
of heavy nuclei (Primakoff effect), e + e _ — > e + e~37r, etc. were suggested for experiments 
0| | ID| . In order then to either deduce the zero energy (the soft pion limit) amplitude from 



experiments or to predict appropriate cross sections from the theoretical side, it should be 
necessary to know how the appropriate physical amplitude may be extrapolated to the zero 
energy limit, or conversely, how the zero energy amplitude be continued to the physical 
region. The original (electromagnetically gauged) Wess-Zumino-Witten term alone does 
not really serve this objective and the standard procedure was to assume the process 
to go via the formation of vector mesons: p (and u) P|pT|. Note that this procedure 
was not unique as will be discussed. Two experimental results have been available: one 
from the Primakoff effect with 40 GeV pion beam ||12|| , see Fig.l, and the other from the 
ix'e -> 7v-7T°e with 300 GeV pions at the CERN SPS @, Fig.2. For both cases, all 
different methods explored in those articles for extrapolation from zero energy amplitude 
fell inside the experimental error bars (about 20% of the central values) once the number 
of colours N c was fixed at three. So one might feel that at least the colour problem 
has found the solution, but that the ATA theorem has not yet been proven: both the 
theoretical models for describing the 7 — > 3n in the zero energy limit and at finite energies 
and the experiments should be improved. In fact, there have been various effective meson 
Lagrangian models constructed so as to be consistent with the underlying QCD (including 
the now fashionable Chiral Perturbation Theory [CHPT], see for example [14]]). And on 
the experimental side (i) there are plans to repeat the pion Primakoff process with 600 
GeV pions at Fermi Lab. (Fig.l), see for example expecting a better statistics and 
making the virtual photons to become closer to the real ones than those in the previous 
experiment [12 1, and (ii) an accepted proposal (at CEBAF) for studying 77r + — > tc + it in 
the photo double pion production: 7^ — > 7i + 7i n (Fig. 3), see |1| using the tagged photon 
beam with energies between 1 and 2 GeV (note that this experiment was investigated 



before [I7|, but the result was strongly in contradiction with the consequences of both 
refs. [0 and [0]). So both theoretical and experimental activities are expected to shed 
a new light on testing the ATA theorem with improved precisions. 

In what follows we shall first discuss some salient features of various effective meson 
Lagrangian approaches for the 77r — > 2tt (or 7 — > 37r) reaction and related processes in 
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Section 2, particularly the ways in which the vector meson contributions are implemented. 
This section may appear a little lengthy, but we consider it useful to give some overview 
of the theoretical trend in a comprehensible manner. Incidentally, by this one might 
see that Chiral Perturbation Theory (CHPT) is not the only relevant approach to low 
energy meson phyiscs. Then in Section 3, we present our confined quark approach for the 
same reaction. This we think should be considered as an alternative to those effective 
Lagrangian approaches discussed in Section 2 where only meson degrees of freedom are 
the relevant quantities (with some minor exception of the appearance of current quark 
mass ratios, quark electric charges, etc.). Section 4 is devoted to an explicit calculation of 
cross section o~/Z 2 for the Primakoff process: the pion production in a nuclear Coulomb 
field on by an incident pion, from some models including our own, and compare with the 



experimental result in |12|]. Then a short conclusion will be given at the end. 



2 Models for Vector Meson Contributions: a Brief 
History 

As mentioned in the previous section, chiral anomaly predicts the zero energy 7 — > 3tt 
amplitude F 3n (0) which is a real constant. It must be continued to appropriate physical 
regions in order to be confronted with relevant quantities extracted from experiments (or 
the continued amplitude must be used to calculate experimentally measurable quanti- 
ties) . The continuation (or extrapolation) then makes the amplitude a function of several 
Lorentz invariant kinematical variables (the number reduces to three, of course, when all 
the particles are physical, of which only two are independent). It was assumed that the 
continuation from zero energy would arise due to the p and to vector meson pole con- 
tributions for reasonable values of these kinematical invariants, say up to ~ 1 GeV 2 Q. 
For the discussion to follow, we choose the process 771 — ► titt to define the kinematical 
variables, see Figs. 1-2. In the physical region for this process, s > 4m 2 and t < should 
hold. Note that in the experiments in [12] |Tj|, the photon is reasonably close to the real 



one: — q 2 < 2 • 10 3 GeV 2 |T2|], and — q 2 m 10 3 GeV 2 [0. For the proposed experiment 



16| , the photon is real: q = but the incident pion with momentum p\ is virtual: p\ < 0. 



Mandelstam variables are defined as 



3 



s = (pi + qf = (P2 +P3) 2 , t = (pi -p 2 ) 2 = (p 3 - q) 2 , u = (pi -p 3 ) 2 = (?2 - <?) 2 , 



so that for processes of [0] and Jl3| s + t + u = "im\ + q 2 . With this preparation we 
shall give a review of several existing models predicting the zero energy (all the particles 
are soft s = t = u = q 2 = 0) amplitude and its continuation to physical regions. As 
stated in the Introduction, ATA ||| found F 3 ^(0) = F 3 ™ from anomalous PCAC 
and current algebra. Then Terent'ev || devised a simplest continuation from the chiral 
anomaly prediction in the form 



F s Js,t,u) 



3tt 



1 + C p e 1 ' 



+ 



+ 



u 



u 



m:. 



(4) 
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3 pane 



0.434, 



nr 



where 5 and 5' are unknown relative phases. This is a form inferred purely intuitively 
from the assumptions of (i) the vector meson dominance (hereafter denoted as VMD), and 
(ii) F 3 n(s, t, u) satisfying a once-subtracted dispersion relation with F^ om the subtraction 
constant at zero energy. 

The second approach to be mentioned here is the work of Rudaz ]TJJ based upon 
the pole model of Gell-Mann et al. [ITS] in which VMD is the guiding principle: all the 
radiative meson processes go via productions of vector mesons. First, Rudaz obtained the 
F irn using the SU(3), p universality and the u — > 7r°7 width and found an agreement with 
the prediction from chiral anomaly. Then he calculated the 7 — > 3ir amplitude and found 
that without anomalous PCAC the ATA theorem could be obtained only if the revised 



Kawarabayashi-Suzuki-Riazuddin-Fiazuddin (KSRF) relation Jl9| would be adopted (a la 
Basdevant and Zinn- Justin (BZ) p0| ) to relate the vector meson and pion quantities at 
low energies: by defining 



Oik = 



q 2 F 2 

mi 



(5) 



with m p and g p7T7T the p meson mass and universal coupling constant (the universality 
means, for example, g p7T7T = q p nn) respectively, the KSRF relation fixes its value as 
«fc = 1/2 (its current empirical value is about 0.53) whereas BZ obtained = 1/3. With 
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this latter value the resultant Rudaz amplitude reads 

F 3 „(s, t, u) = F—^^^i^- + + -^—) ■ (6) 

m t — Q o xrrip — s m p — t m A p — uj 

As a function of s the ratio F 3n (s, 0, 0)/F^ om from the Terent'ev (with unknown phases 

set equal to zero) and the Rudaz models are not identical but the difference is seen as 

at most 10%, which might roughly define what kind of experimental accuracy would be 

required if one wished to distinguish the model difference from given data. 

The story became quite interesting when Rudaz later claimed that in view of the 

later development the above BZ choice: = 1/3, would be somewhat dubious, and that 

the KSRF choice: = 1/2 must be adopted. Then in order to respect the ATA theorem 

he suggested to introduce a non-vanishing direct contact term for the uj — ► 37r coupling 

with its value to be fixed at G u = —g p7T7T /167i 2 F^, in the zeroth order in m 2 /m 2 . This led 

to the revised Rudaz amplitude: 



77? 

F^{s ) t ) u)=F*r " 



2 



1 3 m 2 ,/ 1 J 1 



2 2 3 \m 2 — s m 2 — t m 2 — u 



(7) 



Thus the naive pole contribution in VMD model alone appeared insufficient and an empir- 
ical rather large and negative direct 73tt contact term contribution seemed indispensable. 
This result was later given a support by Cohen [^] who recalled the structure of the 
anomalous Ward identity for 7 — > 3tt discussed by Aviv and Zee from which the 
sign and magnitude of the contact and vector meson pole contributions for the above 
expression could be read off. 

In relation to the forgoing discussion, it may be important here to touch upon the 
subject of the possible non-zero co>3tt contact term the history of which is rather old, but 
still lives in modern effective meson Lagrangian approaches. Earlier, in a pole dominance 
approximation to the isospin one pre scattering combined with the assumption of unsub- 
tracted dispersion relation p3| , its possible existence was suggested once a consistency 
with Bjorken limit was imposed. Later, there was an attempt to determine its magnitude 
from the point of view of hard current algebra plus anomalous Ward identity for three 
and four point functions in describing the u) — ■> 7t + 7t~7t° and 77T — > 2tt ||24j| , but no definite 



value could be reached. Within the experimental accuracy available |[25| , the uj — > 3tt 
transition could be explained by the pole model of VMD |18[ with no apparent need for 
the finite contact term. So for a while the subject stayed dormant. 
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From mid-80's there arose interest in constructing effective meson Lagrangian by 
adding vector (pseudo-vector) meson fields to the combination of the ordinary pseu- 
doscalar (chiral) Lagrangian and the Wess-Zumino-Witten (WZW) term. This was a 
revival of the once successful vector meson dominance model (VMD) in a more refined 
form inferred from underlying QCD. Of those we single out the works by Kaymakcalan et 
al. Pfl, Brihaye et al. ||, and Kuraev and Silagadze Kaymakcalan et al. H] tried a 
non-abelian (vector and axial vector fields) gauging of the chiral plus WZW action and, 
after introducing the VMD (to couple the photon with vector mesons), they required to 
arrive at what is called the Bardeen form of the anomaly. This was found to correctly 
predict the 7r° — > 27 decay rate through the upn coupling: a kind of prerequisite for 
any sensible model. The model has given a definite value for the uj3tt contact term as a 
function of ol\~ as introduced before. For the KSRF choice: = 1/2, which has given 
about 1/4 of the Rudaz value for the contact term, this model did successfully reproduce 
the u) —>■ 3n rate but did not furnish the ATA theorem (in fact, Rudaz has stated JTI]] that 
for any sensible value of the low energy theorem is not reproduciable in this model). 
Note, however, that as an apparent converse, the Rudaz value for the contact term cannot 
supply the sufficient u — > 3tt rate consistent with experiment [36|. It should be interest- 
ing to remark in passing that the authors of ref. |26|] have demonstrated that once the 
kaon loop correction with the WZW term for the K + K~ — > 3tt contact interaction to the 
7 — > 37r amplitude is added within the model of ref. || , the ATA theorem can be restored 
upon suitable regularisations, the consistency condition for which makes = 0.55, very 
close to the empirical value (the pion loop contribution is shown to be small O(m^), thus 
may be disregarded). 



The second set of models @ f26| demanded only the U(l) electromagnetic gauging of 



the combined chiral plus WZW effective action, then tried to consistently introduce vector 
mesons in an empirically well-established principle of VMD. The models easily passed the 
test to reproduce the tt° — > 27 rate. Here again the parameter enters and with its 
value set equal to the KSRF one, the u3tt contact term agrees with the Rudaz one, so 
the ATA theorem was shown to be satisfied (this means, however, that the prediction for 
the u — > 3tt turned smaller than the empirical value). Now summing up, those and other 
effective Lagrangian approaches may need to be confronted with other radiative meson 
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processes, etc. and further relevant experiments are desirable. We note in this regard 
that a discussion in should be quite useful. 



Last in line to be discussed is the chiral perturbation theory (CHPT) applied to 
7T7 — > 2tt and rj — > 27r7 by Bijnens et al. P3. They started with the same chiral plus 
WZW action like in the works discussed above |J |5| [f2B] and, calculated the corrections 
to F^ om which, as has already been stated previously, is obtained upon simply gauging 
the WZW term by the photon field. The corrections include one loop diagrams involving 
one vertex from the WZW term, and tree diagrams from the 0(p 6 ) effective Lagrangian 
appropriate for the process which is required also to cancel the divergence from the loop- 
WZW combination. The finite part of the 0(p 6 ) tree coefficients were fixed by assuming 
their saturation by the nonet vector meson contributions and VMD, then by integrating 
over the vector meson fields (so, needless to say, no explicit vector mesons but only the 
original pseudoscalar octets appear in the effective Lagrangian). The resultant correction 
term is shown to have the following form, 

Fg HPT (s, t, u) = F™ (l + C£js, t, u) + 3 S), (8) 



2w? p J' 



where C^ opg (s,t,u) is the loop correction: 



™° / , x rn 2 n f ml s-A Is -4 l + \/l-A/s 

Ci O0 Js,t,u) = (In J -7- —J— m 



1- -4/5 



t-A /t-4, Vl-4/t + l u-A Fa^4, Jl-A/u+l 
-^ln-^p= = 1 \ —In v / = 

~t Ji-A/t-1 3 V ~ u Jl-A/u-1 



s - A s — 4 
—in \ 

3 V s 

Here s = s/m^, i = t/m%, and u = ujm\. The appearance of the m p in logarithm 
is a standard feature arising from the mass scale needed for dimensional regularisation 
adopted: a common practice in CHPT is to set this scale equal to the p mass. Needless 
to say that the 0(p 6 ) correction should vanish in the chiral limit. It may be of use to 
mention that the origin of the third term in the big bracket in Eq.(8) due to the 0(p 6 ) 
tree contribution is basically the second term in the linearised vector meson propagator 
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contributions, viz. 



— s m p 



together with s + t + u = 3m 2 , (q 2 = for the physical photon). This approximation may 
work typically up to \s,t or u\ w 8m^.: roughly the validity limit of the CHPT expansion. 
For a physical 77r — > 2tt process (where s > 4m 7r 2) this linearised version is therefore 
expected to behave differently from the VMD results discussed above beyond this limit 
for s, even without considering the loop contribution. 

3 A Confined Quark Approach to the 7 — > 3tt Ampli- 
tude 

In the previous section we have seen effective meson Lagrangian approaches which were 
used not only to find the relation between F nll and F 3n (0): to establish the ATA low 
energy theorem, but to continue the latter quantity into appropriate physical regions. 
All modern versions of these approaches are based upon the sum of the lowest order 
chiral Lagrangian and the corresponding Wess-Zumino-Witten (WZW) term in order to 
ensure the possibility of retaining chiral anomaly which allows transitions: even number 
of pseudoscalar mesons <-> odd number of pseudoscalar mesons. The way in which the 
photon and vector (and axial vector) meson fields are introduced to this combination 
differs in different models: some are rather close to each other while others are not. But 
the basic guideline (often employed in implementing the vector-axial vector fields) is that 
once these fields are formally integrated out, the resulting action functional is just the 
original chiral plus WZW term gauged electromagnetically. Thus in a way those effective 
meson Lagrangians are consistent with chiral anomaly by construction and in one way or 
another (recall the kaon loop correction [25| for the model of ref. M) should reproduce the 



ATA theorem. In this regard we have been curious about if there may be a way to approach 
our chirally anomalous processes from the quark model point of view. Our interest in this 
regard is not in those models like in |29j or ]30| in which either specific quark models or 



QCD Lagrangian is adopted in attempting to derive the effective (non-renormalisable) 
meson Lagrangian. Rather, we want to regard mesons as explicit quark- ant iquark bound 
states, so photons couple to mesons through quarks, and various electromagnetic and 
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mesonic transitions proceed explicitly through quark loops. Thus our vector mesons are 
not of Yang-Mills (non-abelian gauge) type. Besides, no VMD assumption is present. 
This may well give different form of amplitudes for the processes of our interest here as 
compared with those from the effective meson Lagrangians discussed so far. In principle, 
any quark model which may be able to describe low energy meson properties will do, 
and here we shall adopt a confined quark model: Quark Confinement Model developed 
at Dubna (QCM). There is a detailed monograph |3lJ on this approach so we only supply 
here the minimum prerequisite for the present purpose. 

In QCM [[H]], generally there are contributions both from quark loops and meson 
pole diagrams. The gauge invariance of the Lagrangian under electromagnetic interaction 
leads to the gauge invariant form of the transtions V — > 7 (V: vector meson) that gives 
a vanishing contribution with real photons. The model is specified by the interaction 
Lagrangian describing the meson transition into quarks, for example, for ir and p mesons 
one can write 



L 7 (x) = ^vf(x)g(x)f 7 5 g(x) + ^=p^(x)q{x)T^q(x) + eA^(x)q(x)Q^q(x) (9) 
v 2 v 2 

with the meson-quark-quark coupling constants gM determined by the compositeness condition: 
that the renormalisation constants of the meson fields are equal to zero so no elementary 
meson: 



'M 



1 - g 2 M W M (m 2 M ) = 0, (10) 



where the empirical meson mass is the input. Here U' M is the derivative of the meson 
mass operator. This way of determining the quark- meson coupling strength is one of the 
differences as compared with other approaches using quark loops like ref.|J2| in which 



actual processes like 7r° — > 2 7 is fit to determine the coupling. Note that the above 
equation provides the right normalization of the charge form factor -Fm(O) = 1. This 
could be readily seen from the Ward identity 

= 9mtLp"^^- = 9m^P^(p, P ) = F M (0) = 1. 
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where T^ I (p,p') is the three-point function describing the meson electromagnetic charge 
form factor. It may be of use to add that the compositeness condition ensures that unlike 
some of the quark-meson hybrid models, we have no danger of double counting. 

Mesonic interactions (meson-photon, meson-meson, etc.) in QCM proceed only through 
closed quark loops: 



da v tr 



M(xi)S v (x 1 - x 2 )...M(x n )S v { 



11^ 



Here. 



^-^/(^"'^ (12) 

is the quark propagator with the scale parameter A characterising the size of the con- 
finement (or equivalent to something like the constituent quark mass), and the measure 
da v , which is essential for quark confinement, is defined to provide the absence of the 
propagator singularities in quark loops, prohibiting the physical quark production: 

" <T " G(z) = a(-z 2 )+zb(-z 2 ). (13) 



v — z 

The shapes of the confinement functions a(u) and b(u), and the scale parameter A 
have been determined from a reasonable model description of the low-energy hadronic 



quantities, see [31] for more details 



a(u) = a exp(— u — a\U) b(u) = b exp(— u +biu). (14) 

Here we have fixed = bo = 2, a\ = l,b\ = 0.4 and A = 460 MeV, which describe 
various basic meson constants quite well, some examples of which may be found in Table 



1, for which various intergals defined in Table 2 are used. See also other quantities in |31 
Note that in Table 2 and in the rest of this work the dimensionless meson masses defined 
as Hm = m M /A appear. We also note that some of the functions in the latter table 
define the meson-quark-antiquark coupling strength through the compositeness condition 
in Eq.(10): for example 



3<£ 2 3a 



p 



47T 2 i?Pp(/4) 4vr 2 Rvv(j% 
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Naturally, the coupling constants g p7r7r , 57^7, and g P7 obtained in the model describe 
the decay widths of the p successfully, which read 



T(p — > 7T7r) 



m 



Q 2 

'48^ 



4m 2 I 3 / 2 



7T 7 ) 



j pn-y 



m 2 



, _ 47ra; 2 2 
e + e ) = — — m p c/ p7 , 



Also it may be useful here to mention that from Table 1 and Eq.(5) our model gives 
a k = 0.53 (by adopting the standard value: m p = 770 MeV), identical to the empirical 
value and quite close to the KSRF value of 1/2. 

For our present objective we first calculate the n° — > 77 amplitude: F^ n . The contri- 
bution comes only from the triangular quark loop (Fig. 4): no vector meson pole. From 
Table 1 one can read off the expression 



a 



7T77 



2R PVV {ti)R P {ti) 



0.94- 



a 



7cF n 



0.024 GeV 



-1 



(15) 



in excellent agreement with the experimental value quoted in the introduction and with 
the chiral anomaly prediction: a/iiF^ = 0.025 GeV -1 (here again (j, M = m M /A is used). 
We note that in the chiral limit: m n = 0, the value in the bracket of Eq.(15) becomes 
0.97 even somewhat closer to the chiral anomaly prediction. Thus we may conclude that 
the model is well calibrated to have successfully passed the first hurdle. 

The next step is to obtain F 3w (s, t, u) for 7 — > 37r which, in the context of the present 
model, is described by two types of diagrams in Fig. 5. We will consider an alternative 
physical process 7r 7 — > nn. The invariant matrix element consists of two parts and reads 



M»[n( Pl ) + 7 (g) - n(p 2 ) + 7r(p 3 )] = e^( Pl ) u (p 2 ) ri (p 3 ) K F 3n (s,t,u), 



(16) 



where 



F 3n (s,t,u) 



F^%s,t,u) + F^ c (s,t,u) 



(17) 
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= ±(^y-L{R D ( s ,t) + R n (t,u) + R n (u,s)}, (18) 



F polc (s,t,u) = — 
m. 



nip — s 



(19) 



"p 

The function denning the box-diagram (corresponding to the first diagram in Fig. 5 
which, however, is not drawn as a proper box) is written as 



R n (s,t) = 2 fd 4 a5(l -J2 a i) -a'(-D A (a 
•* i=i L V 



(20) 



where a! is the derivative of function a(u) in Eq.(14), and D4 is equal to 



D 4 (a) = (otict2 + ct2Ct^ + a^a^m^/K +a 2 o;4s/A + aiast/A . 

The value of R n changes very slowly: less than 1% when s and t are varied within 
and near the physical region for the process under consideration. We may therefore safely 
approximate the box contribution by setting s = t = u = neglecting all the mass and 
momentum dependence with quite a good accuracy. Then one trivially finds -Rd(0, 0) = 2 
leading to 

F " X( °' °' 0) * Zp(0)*Jf* = °- 22F ™ (21) 
The vector meson pole (or resonance) contributions: the second diagram in Fig. 5, are 
defined by F polc (s,t,u) in Eq.(19) which contains the pirn and pir^ vertices which are not 
constants but are explicitly dependent upon the Mandelstam variables due to their quark 
loop structure, 



(s) 



2tt 



R 



VPP 



>/A 5 



Rpp^I) \Rw(fi 2 p V 



(22) 



em p 12 R PVV (s/A 2 ) 



A V 3 ^R PP ^l)R vv ^ p ) • 



Fig. 6 demonstrates their dependence on the virtual p meson mass squared. Apparently 
the variations are not quite small. This may provide a rather interesting effect on which 
we shall discuss later. Note that the p-meson decay widths as presented before are defined 
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by the coupling strengths of Eq.(22) with the s- variable set equal to the square of the 
p-meson mass: g p7T7V = # p7r7r (m^), g^ = g P ^(m 2 p ), and g n = g pl {m 2 p ). As is well known, 
there are relations between these constants (see, for example, pfl) 

9 fry 1/ 9p-mri 9pTT-y 7\9pniT ' F-ir-y-y ^ m po o „ ; (23) 

Z e o7T 

which are obtained based upon an assumption on the VMD and PCAC relation. In our 
model the corresponding relations are found numerically as 

9oy l-0"/$W7r> 9fm~i 9p-n-iv ' 0.84e?77.p— - , (24) 

2 e 87r z l < n 

where our result in Eq.(15) has been taken into account. One can see that the relations 
(23) are reproduced within the commonly accepted accuracy of the VMD assumption 
although our model is not based upon it. In this respect we should stress that with 
empirical values for g p7Tn etc. one obtains the result very similar to ours in Eq.(24). 

Upon combining the box and p-pole contributions, we find for null momenta (s=t=u=0, 
soft pions and photon) 

*W(0) = F 37r (0, 0, 0) = F 3 b r(0, 0, 0) + 6 fr™(° W°) (25) 

m p 

[0.22 + 0.83] = 1.05F^ nom , 



4vr 2 F3 

which is very close to the chiral anomaly prediction (the ATA derivation) or, equivalently, 
to the one from the electromagnetically gauged Wess-Zumino-Witten term as discussed 
in the previous section although we have borrowed nothing from those approaches. 

A particularly noticeable feature of our result is the dominant role played by the p 
meson pole contribution, similar to the VMD approaches. In this respect, it may be useful 
to make here a close comparison of our result with the VMD result. Within the VMD 
models (see, for example) |nj the zero energy 7 — > 3tt amplitude from the vector meson 
pole contribution may be obtained as 

Fl M (0) = Sa k ■ F 3 ™, (26) 
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the derivation of which requires the relations in Eq.(23). A possible presence of the direct 
contact term introduces an additional contribution of the form 



F 3 r(0) = T K) • FST, (27) 

thus 

F 37r (0) = FSST + Fl M = b{a k ) + 3a k ] ■ F~- (28) 

While the form of the p meson pole contribution is identical in any VMD approaches, 
that of the contact contribution differs in different models: for example, Refs. |§ find 
l( a k) — 1 — 3a k , whereas |§ obtained 7(a&) = 1 — 3a k + 3/2a| (we recall that the first 
form works well with the ATA theorem while the second can be made consistent with the 
uo — > 3ir width, both with reasonable values of a k ). 

Regarding our result Eq.(25), it may appear quite natural to identify the box contri- 
bution as that coming from the contact term in the effective meson Lagrangian. Then 
although, as we repeat, our model is not based in any way upon the VMD assumption, 
it would be tempting to assume (although somewhat dubious !) that the relations in 
Eq.(23) may well hold among the coupling strengths evaluated at zero energy, viz. among 
g pnn (0), g pny (0), etc. Then from Eqs. (25-28) this seemingly natural but somewhat naive 
assumption leads to a k = 0.28 (from 3a k = 0.83). We note that since the form of the con- 
tact term is not unique, it is wiser not to use it to extract a k . Clearly, the thus obtained 
value of a k appears quite smaller as compared with both the empirical and the KSRF val- 
ues. Does this imply that our model is not consistent with the Vector Dominance Model 
quantitatively, or qualitatively, or even is there anything wrong with our approach? 

In order to investigate this problem, we want to stress that the key point is the 
fact that, as is clear in Eq.(22) in our quark approach, the couplings to the p meson 
are functions of the virtual mass square of the p, and that , as stated before, the mass 
square dependence is not small (see Fig. 6). As a consequence, the p-meson contribution 
in Eq.(25) evaluated at zero p mass square is considerably smaller than the ordinary value 
obtained in VMD models with fixed coupling strengths (corresponding to those strength 
functions evaluated at the physical p mass). This then leads to the extraction of small a k 
mentioned above as long as we do not queation the appropriateness of such an extraction 
and identification as a k . Now if one insists that somehow a consistency with VMD is 
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mandatory in order our approach to be acceptable, it may be possible to reinterpret our 
result for that objective as follows: Eq.(19) which is the p-pole contribution to the 7 — > 3ir 
amplitude in our model may be decomposed into the pure p meson pole and the rest, the 
latter being analytic at the p meson pole, 

— — r — — + R(s) + (s — > t) + (s — > u). (29) 

m p {m p 2 - s) 

The pure p-pole contribution in the above expression is the one appearing as the p meson 
pole contribution in the VMD models. The analytic part is beyond the VMD approach. 
One may then add the analytic part to the box contribution and regard the sum as 
the contact term introduced phenomenologically in the effective chiral vector meson La- 
grangian approaches stated in the last section. With this rearrangement, Eq.(25) may be 
rewritten as 



T? (C\\ panom ( peon , o 2 #P^ ( m p)gp?r7 ( m p) \ 
v rn p r 3iT ' 

where the first term in the above bracket is the contact term contribution defined as 
F con = ^ OX (0 ; 0, 0) 2[g p7Tn (m 2 p )g p7TJ (m 2 p ) - g pn7T (0)g pnl (0)] 

whereas the second term is the standard p pole contribution. It is easy to show that 
this second term may be identified as 3ak of Eq.(5) when relations in Eq.(23) are strictly 
obeyed. In our model where, instead, Eq.(24) holds, the effective value of at is obtained 

as 

af = ^H)9 Z: H) = QM . (9j^W =0 44 . (32) 

WlpPtyK V rfl p J empirical 

Clearly, this is consistent with the KSFR prediction («& = 0.5) within the accuracy 
of the VMD approach. Then from Eq.(31) the corresponding contact term is equal to 
peon _ i Q5 _ fafff — —0.27, which is consistent with the one in effective Lagrangian 
approaches both in sign and magnitude. Thus our result has been shown not in conflict 
with the VMD models. 

Our amplitude with non-vanishing Mandelstam variables may be conveniently re- 
written in terms of F^ om , 



F 37r (s,t,u) = R 



3tt 



0.22 + C p {s)D p {s) + C p (t)D p (t) + C p {u)D p {u) 



(33) 
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where 



C„(s) = -*<™^J»^> and D ( S \~^A 



P^ ' m 3 panom ' P^ > ™2 e ' 

' n p r 3ir " l p b 

4 Comparison of Models With Data from Existing 
Primakoff Process 

Here we first discuss the value of the zero energy 7 — > 3tt amplitude from various models. 
As already discussed, the chiral anomaly prediction is 9.5 (in units of GeV -3 ). Of those 
models incorporating vector mesons, some have reproduced (or have been taylored to 



reproduce) the chiral anomaly value M |26j while some other M gives a larger value 



about 13.0 due to the smaller magnitude of the u37r contact term. The chiral perturbation 
calculation of ref . [ 28[ to 0(p 6 ) as quoted by Moinester |TJj offers the number 10.7 ± 0.5 



(it should be mentioned here that at the zero energy [or the soft pion limit] this approach 
should reproduce the anomaly prediction so the quoted value should be taken with some 
caution). Our QCM yields 10.0 for its value. These values should be compared with the 
one extracted from the Primakoff process |12| is 12.9 ± 0.9 (stat) ± 0.5 (sys). Based upon 
this experimental value it is hard to rule out (or give preference to) any of the prediction. 

Next, the model amplitude will be used to predict the integrated Primakoff cross 
section for the pion production by incoming pions as in Fig.l. We adopt the Mandelstam 
variables as introduced previously and go to the centre of mass system for the two final 
state pions. Then for 77r — > tttt (see, Fig.l) noting that s + t + u = 3m^ + q 2 , we have 

Pi + q = P2 +P3 = 0, cos 9 = - - , Pi = E, 

\Pl\\P2\ 

and taking into account that the photon is almost real (q 2 ~ 0), one may write 
t = ~[3m 2 — s + (s — m^)\j\ — Am 2 /scosd], 

u = -[3m 2 — s — (s — m 2 )^l — Arri^/s cos 9]. 

The physical region for s is Am 2 < s < 16m 2 . We obtain various model amplitudes for 
F 3n (s,t,u) in units of the chiral anomaly prediction: _P 3 a ^ om = e/47r 2 F^ as a function of 
s for the forward scattering angle: Fig. 7 (the angular dependence has been found rather 
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weak so we only give only one angular value). The plotted curves are the original Rudaz 
VMD model result JTTJ], then his modified (improved) amplitude |[21|| , the amplitude from 
Kaymakcalan et al.'s model M, chiral perturbation (CHPT) prediction ||28|| , and our 
present quark model result. Except for the amplitudes from || and other curves 
are mutually rather close. Kaymakcalan et a/.'s || result is due to its smaller sized of 
the contact term which has given a large zero energy amplitude. On the other hand the 
CHPT amplitude above s > 8m 2 has come out quite differently from the others. The 
validity for the CHPT expansion is up to about this value of s, so the resultant difference 
may be natural. If future experiments may be able to extract this s and scattering angle 
dependent behaviour of the amplitude, it would be of quite help. At present, however, 
we should be satisfied with their prediction on the available data [12 . 



The cross section for pion production by the pion incident in the nuclear Coulomb 
field (the Primakoff process) 

7T° + (Z, A) (34) 
• 2tt cross section through the equivalent 



7T 



71 



as in Fig.l may be expressed in terms of the •yn 
photon method: 



where 



and 



da 



Z 2 a 



W 2 



dsdtdq 2 



71 



q 2 - 



q 2 - 

^mm 



da, 



7"7T — >7T7r 



mi 



dt 



s — ml 



2E 



(35) 



da, 



77T — >7T7V 



\F 3 Js,t,u) 



[s — Am 2 ) sin 2 6 



(36) 



dt 512tt 

with the vertex F 37T (s, t, u) defined by the dynamics of the 7 — > 3n transtion. 

Neglecting the g 2 -dependence in F 37r (s, t, u), viz. setting q 2 w (thus s + t + u = 3m 2 ), 
the integrated cross section per unit proton charge in the nucleus is given as 



a 
Z 2 



a 



1024tt 2 



ds [s 4m 2 ) 3 / 2 | ln qL. _ 1 + 0Ln 



invi 



q 2 . 

1 mm 



dd sin 3 6 l |F3 7r (s,t,M) 



(37) 
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The condition for the Primakoff effect is the following Q 







E 2 » ml ~ s » q 2 . 
This condition is better satisfied when the energy of the incident pion beam becomes 



higher. The Serpukhov experiment |12| has been carried out with the 40 GeV pion beam. 
The maximum momentum transfer was g^ ax = 2 • 1(T 3 GeV 2 and s max = 10m 2 . 

The result of the calculation of the cross section for various approaches are shown in 
Table 3. One can see that the results are rather close to each other, except for the one 
from the model of ref. ||, which could be understood because its prediction for F 3n (0) 
is about 30% larger than the chiral anomaly prediction, as discussed in Section 3. We 
note that a somewhat larger value (about 10%) of the cross section in our QCM model 



as compared with those from those of refs. [TJI]] and is also largely the reflection 
that our ^3^(0) is about 5% larger than the chiral anomaly prediction. In Table 3 the 
row denoted as Chiral Anomaly has been obtained from replacing F 3n (s,t,u) simply by 
F 3 ™ m = e/4n 2 F^, thus disregarding the energy dependence (and angular variations) in 
calculating a/Z 2 above. This is found to reduce the integrated cross section by more than 
10%. But at present, the extracted values from the experiment may only be able to say 
that the result from ref. || may be inadequate (note, however, that as we have stated 
before, one K-meson loop correction seems to have restored the value of the contact term 
to the Rudaz one |n[]). Improved experiments are definitely in need. 

5 Discussion and Conclusion 

It may sound repetitive, but we make a pedagogical review at the beginning of the dis- 
cussion. 

Only the anomalous (but not the ordinary) PCAC combined with current algebra 
was able to correctly predict the observed 7r° — > 27 decay rate and the zero energy 
7 — > 37r amplitude F 3n (0) was predicted from the former amplitude: we have termed this 
as the ATA theorem || [|| || . This theorem was rederived quite naturally later from the 
electromagnetically gauged Wess-Zumino-Witten action (WZW) || @. Modern versions 
of the vector dominance models (VMD) [|J || |21] |26j were used in an attempt to reach 
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this theorem with the help of the KSRF relation |19| (and its more relaxed form). In 
order to be consistent with the ATA theorem the vector dominance models have been 
required to introduce a non-zero constant four-point coupling for u> — > 3n leading to a 
contact 7 — > 3ir coupling which has an opposite sign as compared with the vector meson 
pole contribution. While the vector meson pole contribution in these models for the zero 
energy 7 — > 3tt amplitude is 1.5 times the value obtained from the ATA theorem, the value 



of the contact term contribution differs: - 1/2 times the ATA value from refs. || [21] [26 
and about 4 times smaller in ref.|| (when evaluated at the tree level), and the former 
is only consistent with the ATA theorem. We have calculated above the Primakoff cross 
section for these models and found that improved experiments may be able to say which 
group of the VMD model (regarding the contact term value) should be more relevant to 
the 7 — > 37r amplitude. 

We have devised a confined quark model (QCM) to study the same subject without 
implementing any explicit chiral anomaly. The result has come out quite consistent with 
the ATA theorem and the consequence of the electromagnetically gauged WZW action. 
The obtained amplitude F 37T (s, t,u) is fairly close to the VMD prediction of refs. [[J ET| 



26fl. However, if we regard our quark box diagram as the contact term contribution, it 



is found to have the same sign as the vector meson pole contribution, and counts just 
22% of the ATA theorem (or the chiral anomaly) prediction: F^ om = e/4n 2 F^. This 
then implies (as it has eventually turned out) that our vector meson pole contribution 
accounts for about 80% of the chiral anomaly prediction. This should be contrasted to 
the result of the VMD prediction discussed above which provides 1.5F^ om upon emptying 
the KSRF relation: = 1/2. At this point we want to remark that, as in Fig. 6 both 
our pnir and p77r coupling vertices are functions of the p meson mass squared arising 
from the quark loop, and that their variations are not negligible, especially the former. 
Their values at m? p are practically the same as those in the VMD models, see Table 1, 
but at zero mass they are smaller. Those smaller values have been used to find the value 
of ~ 80% for the p-pole contribution mentioned above. When those coupling constants 
are evaluated at m 2 p , one finds that the vector meson contribution to F 3n (0) becomes 1.3 
times the chiral anomaly prediction, not very different from the VMD result of 1.5 times 
(recall Eqs.(30,32)). 
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As for the difference in sign and magnitude of the contact term between our quark 
model prediction and those from the VMD models, it does not appear easy to find the 
answer to this discrepancy as long as we stick to regarding the the contact term as our 
box contribution: within our confined quark model approach the loop integral for the box 
diagram is completely convergent and gauge invariant. In the approach in ref.0, the loop 
integral generates a constant term breaking the gauge invariance so it must be subtracted 
out (it is related to the so-called surface term arising from the shifting of the integration 
variable in a divergent integral, see for example section 6.2 of ref.|34|]). This subtraction 
appears to be the origin of the sign change which we do not have in QCM. So what could 
be the possible compromise ? There may be some effects coming from the excited vector 
meson(s) like p' which was expected to cure the trouble with the u — > 3tt rate with the 
larger contact term: — 0.5F 3 a ° om ||27|| . Or the effect of the Regge trajectory contributions 



might be included [26]. Either of these effects might be included in the vector meson pole 
contributions. These lines of considerations might work as a bridge between our and the 
VMD pictures: one might suspect that some type of hadron-quark duality (for example, 
quark loops vs. meson loops) be operative somewhere. 

Once we adopt a somewhat more relaxed definition of our contact term, however, the 
resolution of the apparent discrepancy in sign and the size of the contact term in our and 
VMD models may seem feasible. This could proceed along the line we exploited in Section 
3 when extracting from our prediction of F 37r (0). Namely, we decompose our p-pole 
contribution into the pure pole and analytic parts. Then the analytic part is added to the 
box contribution which one may regard as the effective contact term. Now, this effective 
contact term is no longer a mere constant, but a function of the Mandelstam variables: 
s, t and u. It is then easy to see, as found in Section 3, that in the zero energy limit 
this term obtains a negative sign with magnitude about half the Rudaz value [0 making 
our F^ n (0) prediction very close to the chiral anomaly one . At finite energies (finite s) 
one may be able to see that the magnitude of this term decreases, and this is consistent 
with the apparent need for small (or zero) contect term contribution in the u — > 3tt width 
within the context of the effective chiral vector meson models. We note in this respect 
that both experimental and model investigations, for example, on e + e~ —* 3tt (including 
uj — > 37r) should be quite helpful. 
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From the point of view of the proposed 77r — > tttt experiments (either the photon or 



pion in the initial state is virtual, or at best close to real) ||T5[ |16| to extract -^3^(0) or more 
ambitiously to deduce F 3n (s,t,u) for different s and cm. scattering angle (for example, 
to clarify the distinction between our quark model and VMD predictions: notably the role 
played by either the constant or non-constant contact term, the latter arising from the 
Mandelstam variable dependent p couplings), they will be required to attain a very high 
precision in data acquisition and analyses. In this respect it may be useful to be reminded 
that if some sort of analytic extrapolation to a pole (either for the pion or photon) is needed 
in analysing the data, one must not forget the fact that assumed analyticity property (of 
the relevant amplitude) and the statistical nature of the data are mutually exclusive. So 
instead of a naive straightforward extrapolation, some optimal mapping may have to be 
sought beforehand in order to minimise the unwanted (and in most cases the dominant) 



background contributions, see for example refs.[35| |36| concerning this procedure. Finally, 
in view of the closeness of the predictions from some VMD models, based upon the chiral 
plus the WZW action which basically implements the chiral anomaly, and our confined 
quark model approach (which, at least, is not explicitly built upon chiral anomaly like 
the WZW term), may it be of some use to ask ourselves if the proposed experiments will 
be really testing the ATA theorem as a consequence of chiral anomaly ? 
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Fig. 5. The 7 — > 3n in QCM. The first diagram is called the Box contribution while the 
second type is termed as the vector meson pole contribution. 
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Fig. 7. The 7 — > 3tt amplitude for the physical process 77? — > 27 for zero scattering angle, 
normalised to the chiral anomaly prediction: F^ om = e/AirF^. Various model predictions 
drawn are (1) the VMD model of Rudaz [ TTJ , (2) the modified Rudaz model implementing 



the KSRF relation |21[| , (3) chiral Lagrangian with Wess-Zumino-Witten term gauged by 
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Table 1. Several low energies quantities obtained in QCM. 





QCM 


Expt. 


p _ Ai2 P (/4) 


3 

#pp(/4) 


93 MeV 


93 MeV 


p _ ARp(ji 2 K ) 




112 MeV 


111 MeV 


T? _ e 2 1 RpvviA) 
^ - A v/3^p(^) 


0.024 GeV 1 


0.025 GeV" 1 


2itV2Rvpp(v 2 p ) 


6.0 


6.0 


errip fj R P vv(^ 2 p ) 


0.16 


0.17 






0.18 


0.20 



Table 2. Structural integrals. 



Rpp(x) 


= B + fjdub(-f)^ 


Rw{x) = 


B + ffdub( ?) (1 -^ 2/4) 


R P {x) = 


= A + f fdua(-fWl-u 


Rpvvix) 




Rvpp{x) 


= S + f /d«6(-f)>/l-« 


Rv(x) = Bo 


+ f/«(-^)(l + w/2)v / l-w 




oo 

Aq = J dua(u) 


= 1.09 B = Jdub{u) = 2.26 
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Table 3. 



a/Z 2 (nb) 


Model 


Expt. 12 


Chiral anomaly 


0.95 




Modified Rudaz [|2Tl 


1.17 


1.78 ±0.47 (C) 


Bijnens et al. 28 


1.12 


1.54 ±0.34 (Al) 


Kaymacalan et al. [|8j 


2.10 


1.64 ±0.37 (Fe) 


QCM 


1.27 
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